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ABSTRACT 


Two-dimensional  continuum  modeling  and  simulations  were  conducted  to  predict 
how  the  size,  quantity,  and  stiffness  of  reinforcing  particles  such  as  carbon  nanotubes 
(CNTs)  affect  failure  mechanisms  at  the  interface  of  composite  structures.  First,  the 
strength  model  used  the  finite  element  method  (FEM)  on  a  slender  composite  beam  with 
step-joint  containing  reinforcing  particles  to  predict  its  critical  stress-strain  behavior  at 
the  joint  interface  under  compressive  axial  load.  Next,  the  fracture  mechanics  model 
used  the  virtual  crack  extension  method  on  the  same  composite  beam  containing  an 
internal  crack  to  predict  how  the  energy  release  rate  was  affected  by  reinforcing  particles 
at  the  interface  under  the  same  compressive  axial  load.  Comparing  the  two  results  to 
experimental  data  showed  that  the  fracture  mechanics  model  predicted  the  interface 
failure  behavior  better  than  the  strength  model.  Finally,  the  fracture  mechanics  model 
was  used  for  a  composite  plate  containing  an  edge  crack  to  study  how  the  energy  release 
rate  was  affected  by  several  parameters  of  reinforcing  particles  near  the  crack  tip  under 
transverse  shear  load.  In  each  case,  homogeneous  models  served  as  baselines  for 
comparative  analyses.  Outcome  of  this  work  not  only  represents  reliable  and  efficient 
modeling  of  composite  interfaces  in  order  to  improve  failure  strength  through  the 
addition  of  nanoscale  reinforcing  particles  such  as  CNTs  but  also  serves  to  focus  future 
research  in  structural  application  of  CNTs,  especially  within  testing  and  evaluation  of 
CNTs  in  composite  scarf-joint  interfaces. 


v 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


vi 


TABLE  OF  CONTENTS 


I.  INTRODUCTION . 1 

A.  BACKGROUND . 1 

B.  OBJECTIVES . 2 

II.  CONTINUUM  MODELING . 3 

A.  INTRODUCTION . 3 

B.  STRENGTH  MODEL . 3 

1.  Cantilevered  Beam . 3 

2.  Critical  Load  Beam  Bending . 4 

3.  FEM  Model  of  Composite  Beam . 4 

C.  FRACTURE  MECHANICS  BEAM  MODEL . 5 

1.  Cantilevered  Beam  with  Internal  Crack . 5 

2.  FEM  Model  of  Composite  Beam  with  Internal  Crack . 6 

3.  Virtual  Crack  Extension  Method . 7 

4.  Energy  Release  Rate . 9 

D.  FRACTURE  MECHANICS  PLATE  MODEL . 9 

1.  Composite  Plate  with  Surface  Crack . 9 

2.  FEM  Model  of  Composite  Plate  with  Surface  Crack . 10 

3.  Virtual  Crack  Extension  Method . 1 1 

4.  Energy  Release  Rate . 13 

III.  RESULTS  AND  DISCUSSION . 15 

A.  STRENGTH  MODEL . 15 

1.  Homogeneous  Beam . 15 

2.  Baseline  Composite  Beam  without  Reinforcing  Particles . 15 

3.  Composite  Beam  with  Reinforcing  Particles  in  Matrix  Material..l6 

B.  FRACTURE  MECHANICS  BEAM  MODEL . 20 

1.  Baseline  Composite  Beam  without  Reinforcing  Particles . 20 

2.  Composite  Beam  with  Reinforcing  Particles  in  Matrix  Material..21 

C.  FRACTURE  MECHANICS  PLATE  MODEL . 22 

1.  Homogeneous  Plate . 22 

2.  Baseline  Composite  Plate  without  Reinforcing  Particles . 22 

3.  Composite  Plate  with  Reinforcing  Particles  in  Matrix  Material. ..23 

4.  Composite  Plate  with  Varying  Particle  Position . 23 

a.  Particle  at  position  a . 24 

b.  Particle  at  position  b . 24 

c.  Particle  at  position  c . 25 

d.  Particles  at  positions  a,  b,  c . 26 

5.  Composite  Plate  with  Varying  Particle  b  Distance  to  Crack  Tip  ..29 

6.  Composite  Plate  with  Varying  Particle  Size  at  Prescribed 

Distance  to  Crack  Tip . 31 

7.  Composite  Plate  with  Varying  Particle  Stiffness  Nearest  Crack 

Tip . 33 

vii 


IV.  CONCLUSIONS  AND  RECOMMENDATIONS . 37 

A.  CONCLUSIONS . 37 

B.  RECOMMENDATIONS . 38 

LIST  OF  REFERENCES . 39 

INITIAL  DISTRIBUTION  LIST . 41 


viii 


LIST  OF  FIGURES 


Figure  1 .  Composite  beam  under  applied  compressive  load . 4 

Figure  2.  2-D  FEM  8-node  structural  solid  element  geometry  (PLANE183)  [14] . 4 

Figure  3.  Strength  model  of  composite  beam  under  compressive  nodal  load . 5 

Figure  4.  Composite  beam  with  internal  crack  under  applied  compressive  load . 6 

Figure  5.  Mesh  of  fracture  mechanics  beam  model  with  internal  crack  ( a )  under 

compressive  nodal  force  with  expanded  view  of  crack  tip  (b) . 7 

Figure  6.  Schematic  interpolation  of  crack  tip  forces  &  displacements  for  fracture 

mechanics  beam  model . 8 

Figure  7.  Composite  plate  with  surface  crack  under  applied  shear  load . 10 

Figure  8.  Mesh  of  fracture  mechanics  plate  model  with  surface  crack  (a)  under  shear 

nodal  force  with  expanded  view  of  crack  tip  ( b ) . 11 

Figure  9.  Schematic  interpolation  of  crack  tip  forces  &  displacements  for  fracture 

mechanics  plate  model . 13 

Figure  10.  Composite  beam  matrix  material  containing  reinforcing  particles  of  three 
sizes  A  &  two  volume  fractions  ® :  A0=7%  =  1 .23 mm2  (a), 

=  1  ■23/h/h2  (b),  Alll=7%  =  4.91mm2  (c),  A0=l4%  =  4.91mm2  (d), 

A>=n  =  1 1  -0mm2  (e),  A0=14%  =  1 1.0mm2  (f) . 17 

Figure  1 1 .  Combined  maximum  von  Mises  equivalent  stresses  cre  in  matrix  material 

of  composite  beam  with  reinforcing  particles  for  volume  fractions 

®  =  7%  &  ®  =  14% . 19 

Figure  12.  Average  maximum  von  Mises  equivalent  stress  <Je  in  matrix  material  of 
composite  beam  with  reinforcing  particles  at  volume  fraction  ®  =  7% , 

®  =  14% ;  comparison  to  base  model . 20 

Figure  13.  Composite  plate  matrix  material  containing  reinforcing  particles  a,  b,  c  in 

position  radial  distance  r  to  crack  tip . 23 

Figure  14.  Combined  energy  release  rates  Gn  for  reinforcing  particles  a,  b,  c,  in 
matrix  material  of  composite  plate  varying  particle  stiffness  from 
E3=8.0GPa  to  E3  =500GPa  with  particle  size  A  =  0.03961mm2  at 

distance  r  =  1.124mm  to  crack  tip . 28 

Figure  15.  Average  energy  release  rate  Gn  for  reinforcing  particles  a,  b,  c,  in  matrix 
material  varying  particle  stiffness  from  E3  =  tt.QGPa  to  E3  =  5()()GPa 

with  particle  size  A  =  0.0396 Imm2  at  distance  r  =  1.124mm  to  crack  tip; 

comparison  to  baseline  model . 29 

Figure  16.  Energy  release  rate  Gn  for  reinforcing  particle  b  in  matrix  material  of 
composite  plate  varying  distance  r  to  crack  tip  from  r  =  0  to  r  =  3 mm 

with  particle  stiffness  E3  =  400 GPa . 31 

Figure  17.  Energy  release  rate  Gn  for  reinforcing  particle  b  in  matrix  material  of 
composite  plate  varying  particle  size  from  A  =  0.01mm2  to  A  =  2mm1 


IX 


with  particle  stiffness  E3  =400GPa  at  prescribed  distance  r  =  1.124mm 

to  crack  tip . 33 

Figure  18.  Energy  release  rate  Gn  for  reinforcing  particle  b  in  matrix  material  of 
composite  plate  varying  particle  stiffness  from  E3=8.0GPa  to 

E3  =  500 GPa  with  particle  located  nearest  the  crack  tip  at  r  =  0.3mm . 35 


x 


LIST  OF  TABLES 


Table  1 .  Maximum  stress  and  strain  of  homogeneous  beam  ( E  =  50 GPa  ) . 15 

Table  2.  Maximum  stress  and  strain  in  matrix  material  of  baseline  composite  beam 

where  Ex  =  50 GPa  ,  E2  =  8.0 GPa  ,  without  reinforcing  particles . 16 

Table  3.  Maximum  stress  and  strain  in  matrix  material  of  composite  beam  where 
E]  =  50GPa  ,  E2  =  S.OGPa  ,  E3  =  400 GPa  ,  varying  particle  size  from 

A  =  1.23mm2  to  A  =  11.0mm2  while  holding  constant  volume  fraction  ® 

at  either  of  ®  =  7%  or  ®  =  14% . 18 

Table  4.  Average  increase  in  maximum  vonMises  equivalent  stress  ae  varying 
reinforcing  particle  size  A  in  matrix  material  of  composite  beam  when 

compared  to  baseline  model . 18 

Table  5.  Energy  release  rate  Gu  in  matrix  material  of  baseline  composite  beam 

where  Ex  =  50 GPa  ,  E2  =  8.0 GPa . 21 

Table  6.  Energy  release  rate  Gu  in  matrix  material  of  composite  beam  with 

reinforcing  particles  where  E2  =  400 GPa . 21 

Table  7.  Energy  release  rate  Gu  in  homogeneous  plate  ( E  =  50 GPa  ) . 22 

Table  8.  Energy  release  rate  Gu  in  matrix  material  of  baseline  composite  plate 

where  Ex  =  50 GPa  ,  E2  =  8.0 GPa  without  reinforcing  particles . 22 

Table  9.  Energy  release  rate  Gu  for  reinforcing  particle  a  in  matrix  material  of 
composite  plate  varying  particle  stiffness  from  E:,  =  8.0  GPa  to 

E2  =  500 GPa  with  particle  size  A  =  0.03961 mm2  at  distance  r  =  1.124 mm 

to  crack  tip . 24 

Table  10.  Energy  release  rate  Gu  for  reinforcing  particle  b  in  matrix  material  of 
composite  plate  varying  particle  stiffness  from  E:,  =  tt.QGPa  to 

E:.  =  500 GPa  with  particle  size  A  =  0.03961 mm2  at  distance  r  =  1 . 1 24 nun 

to  crack  tip . 25 

Table  11.  Energy  release  rate  Gu  for  reinforcing  particle  c  in  matrix  material  of 
composite  plate  varying  particle  stiffness  from  E:,  =  ti.QGPa  to 

E:.  =  500 GPa  with  particle  size  A  =  0.03961 mm2  at  distance  r  =  1.124 mm 

to  crack  tip . 26 

Table  12.  Energy  release  rate  Gu  for  reinforcing  particle  a  be  in  matrix  material  of 
composite  plate  varying  particle  stiffness  from  E:,  =  tt.QGPa  to 

E2  =  500 GPa  with  particle  size  A  =  0.03961 mm2  at  distance  r  =  1 . 1 24mm 

to  crack  tip . 27 

Table  13.  Average  decrease  in  energy  release  rate  Gu  for  reinforcing  particles  a,b,c, 

when  compared  to  baseline  model . 27 


xi 


Table  14.  Energy  release  rate  Gn  for  reinforcing  particle  b  in  matrix  material  of 
composite  plate  varying  distance  r  to  crack  tip  from  r  =  0  to  r  =  2.937 mm 
with  particle  stiffness  E3  =  400 GPa . 30 

Table  15.  Average  decrease  in  energy  release  rate  Gn  for  particle  b  varying  distance 

r  to  crack  tip  when  compared  to  baseline  model . 30 

Table  16.  Energy  release  rate  Gn  for  reinforcing  particle  b  in  matrix  material  of 
composite  plate  varying  particle  size  from  A  =  0.01mm2  to  A  =  2mm2 
with  particle  stiffness  E:,  =  400 GPa  at  prescribed  distance  r  =  1.124mm 

to  crack  tip . 32 

Table  17.  Average  decrease  in  energy  release  rate  Gn  for  particle  b  while  varying 

size  A  when  compared  to  baseline  model . 32 

Table  18.  Energy  release  rate  Gn  for  reinforcing  particle  b  in  matrix  material  of 
composite  plate  varying  particle  stiffness  from  E^=8.0GPa  to 

E 3  =  500 GPa  with  particle  located  nearest  the  crack  tip  at  distance 

r  =  0.3  mm . 34 

Table  19.  Average  decrease  in  energy  release  rate  Gn  for  particle  b  while  varying 

stiffness  when  compared  to  baseline  model . 34 


xii 


I.  INTRODUCTION 


A.  BACKGROUND 

Since  first  being  reported  by  Iijima  in  1991  [1],  carbon  nano  tubes  (CNTs)  have 
been  the  subject  of  widening  research  in  fracture  mechanics  given  their  exceptional 
mechanical  properties.  With  elastic  modulus  on  the  order  of  lTPa,  CNTs  have  potential 
for  a  multitude  of  structural  applications  including  use  as  reinforcing  particles  in  polymer 
matrix  composites  or  perhaps  even  metal  matrix  composites  in  order  to  enhance  the  joint- 
interface  fracture  toughness  [2-3]. 

A  survey  of  recent  work  showed  classical  molecular  dynamics  (MD)  simulations 
had  been  used  to  model  elastic  and  plastic  properties  of  CNTs  as  well  as  mechanisms  of 
strain  release  under  various  types  of  loading  conditions.  In  one  case  [4],  MD  simulations 
found  CNTs  to  be  exceptionally  resilient  and  able  to  sustain  extreme  strain  deformation 
with  no  signs  of  brittleness,  plasticity,  or  atomic  rearrangements.  In  another  case  [5], 
both  MD  simulations  and  experimental  testing  were  combined  in  continuum  modeling 
that  revealed  single-  and  multi-walled  CNTs  are  highly  flexible,  able  to  withstand  large 
strains  and  even  bending  in  excess  of  110°.  In  fact,  ref.  [6]  used  MD  simulations  to 
perform  a  comparative  analysis  of  elastic  modulus  on  two  common  types  of  nanotubes, 
suggesting  there  exists  a  58%  greater  stiffness  in  single-walled  when  compared  to 
bamboo  structure  CNTs, 

Remarkable  achievements  have  also  been  made  in  the  application  of  CNTs  to 

structures,  most  notably  through  CNT-reinforced  polymer  matrix  composites  [7-9]. 

Uniform  dispersion  within  the  polymer  matrix  and  wettability  remain  critical  issues.  In 

one  experiment  [10],  the  addition  of  1  wt.  %  CNTs  to  polystyrene  resulted  in  a  favorable 

42%  increase  in  stiffness  as  well  as  a  25%  increase  in  strength.  In  another  experiment 

[11],  a  single  multi- walled  nano  tube  was  embedded  in  a  polymer  matrix  was  found  to 

have  very  high  compressive  strength,  more  than  2  orders  of  magnitude  higher  than  any 

known  fiber.  In  fact,  ref.  [12]  dispersed  7.5  g/m  of  various  multi-walled  and  bamboo 

CNTs  into  a  composite  scarf  joints  containing  vinyl  ester  matrix  and  subjected  them  to 
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compressive  axial  loading  that  resulted  in  a  26%  increase  in  stiffness  as  well  as  an  11% 
increase  in  strength  on  average  when  compared  to  similar  joints  without  CNTs. 

But  while  design  and  fabrication  of  CNTs  is  clearly  on  the  nanoscale,  there  exists 
a  fundamental  need  for  favorable  macroscopic  results  in  engineering  applications.  It  is 
through  the  use  of  continuum  modeling  that  the  seemingly  theoretical  design  and 
fabrication  of  nanoscale  particles,  such  as  CNTs,  can  be  readily  applied  to  engineering 
designs  in  order  to  predict  macroscopic  mechanical  behavior  of  a  system.  This  type  of 
modeling  is  essential  to  the  development  of  CNT-reinforced  composites. 

B.  OBJECTIVES 

The  principal  objective  of  this  work  is  to  provide  a  reliable  and  efficient  method 
of  modeling  composite  interfaces  in  order  to  improve  failure  strength  through  the 
addition  of  nanometer  scale  reinforcing  particles  such  as  CNTs.  To  this  end,  modeling 
will  be  perfonned  in  order  to  predict  how  the  size,  quantity,  and  stiffness  of  reinforcing 
particles  such  as  CNTs  affect  failure  mechanisms  at  the  interfaces  of  composite 
structures.  Comparisons  will  be  drawn  to  recent  experimental  work  in  e.g.,  [12].  The 
primary  goal  here  is  to  assist  future  researchers  in  any  structural  application  of 
nanotechnology,  especially  within  the  realm  of  testing  and  evaluation  of  CNTs  in 
composite  scarf-joint  interfaces. 
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II.  CONTINUUM  MODELING 


A.  INTRODUCTION 

Two-dimensional  modeling  can  be  used  to  predict  the  fracture  mechanics  of 
composite  material  reinforced  by  particles  such  as  carbon  nanotubes  (CNTs).  Varying 
the  size,  location,  and  stiffness  of  reinforcing  particles  at  the  interface  causes  changes  in 
stress-strain  behavior  as  well  as  the  amount  of  energy  released.  Identifying  critical  values 
for  each  helps  in  the  prediction  of  crack  propagation  and  provides  a  correlation  between 
the  analytical  results  presented  here  and,  for  example,  the  experimental  results  obtained 
in  [12]. 

B.  STRENGTH  MODEL 

The  strength  model  shows  how  the  stress-strain  behavior  is  influenced  by  the 
application  of  reinforcing  particles  in  the  matrix  material  of  a  composite  beam.  These 
analytical  predictions  are  made  by  calculating  the  state  of  stress  and  strain,  and  related 
energies,  at  the  interface  in  order  to  predict  failure  under  an  applied  compressive  load. 

1.  Cantilevered  Beam 

Consider  a  composite  beam  as  shown  in  Figure  1.  This  slender  cantilevered  beam 
comprises  three  linear  isotropic  materials  having  Young’s  moduli  Ex=50GPa, 

E2  =  8.0 GPa  ,  E3  =  400GPa  and  Poisson’s  ratios  v,  =  0.30  ,  v2  =  0.28  ,  v3  =  0.20  .  Let 
length  /  =  0.62 m  and  height  h  =  0.05 m  .  The  left  edge  is  clamped  while  the  right  edge  is 
subject  to  compressive  force  Fx  =  5,000V  and  constrained  in  the  y-direction.  To  control 

for  the  effects  of  bending,  surface  nodes  at  the  loading  point  are  coupled  in  the  y- 
direction  thus  preventing  the  free  end  from  rotating  while  still  allowing  for  axial 
displacement. 
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Figure  1. 


Composite  beam  under  applied  compressive  load. 


2.  Critical  Load  Beam  Bending 


Using  first  principles  of  Euler  elastic  stability  for  slender  beam  to  analyze  the 
critical  load  Pcr  of  the  main  constituent  material 


=  203.1-106Ar  (1) 


where  the  effective  length  Le  =  0.70  (/)  for  clamped  end,  results  in  Pcr  =  203 .1-1  O'1  N  , 

which  is  substantially  greater  than  the  load  P  =  5kN  being  applied  here  thus  buckling  is 
not  expected  [13]. 

3.  FEM  Model  of  Composite  Beam 

A  two-dimensional  finite  element  method  (FEM)  model  in  ANSYS  uses  higher 
order  eight-node  element  (i.e.,  PLANE  183)  having  quadratic  displacement  behavior. 
Each  node  has  two  degrees  of  freedom,  ux  and  uy,  with  the  degenerated  element  having 
coincident  nodes  K,  L,  and  O,  as  depicted  by  the  shape  in  Figure  2  [14].  Linear  elastic 
plane  stress  behavior  is  assumed  for  this  static  analysis. 


v 

(<x  axial) 


X(or  radial) 


Figure  2.  2-D  FEM  8-node  structural  solid  element  geometry  (PLANE183)  [14]. 


Each  area  is  freely  meshed  using  smart  size  level  1  elements  with  one  subsequent 
refinement.  The  right  side  of  the  beam  has  an  applied  total  force  load  Fx  =  -5kN 
distributed  to  each  node  and  constrained  in  the  y  direction.  The  left  side  of  the  beam  is 
clamped  and  therefore  constrained  in  both  x  andy  directions.  See  Figure  3. 


Figure  3.  Strength  model  of  composite  beam  under  compressive  nodal  load. 

C.  FRACTURE  MECHANICS  BEAM  MODEL 

The  fracture  mechanics  beam  model  shows  how  an  internal  crack  or  flaw  in  the 
interface  of  a  composite  joint  propagates  under  an  applied  axial  load  and  how  this  process 
is  influenced  by  the  application  of  reinforcing  particles  such  as  CNTs.  These  analytical 
predictions  are  made  by  calculating  fracture  parameters  such  as  the  energy  release  rate  in 
the  region  of  a  crack  in  order  to  estimate  the  crack  growth  rate.  While  the  energy  release 
rate  G  represents  the  amount  of  work  associated  with  a  crack  opening,  it  is  at  the  point  of 
instability  when  the  change  in  potential  energy  equals  the  change  in  the  crack  length  that 
fracture  occurs. 

1.  Cantilevered  Beam  with  Internal  Crack 

Consider  a  composite  beam  with  internal  crack  as  shown  in  Figure  4.  This 
slender  cantilevered  beam  comprises  three  linear  isotropic  materials  having  Young’s 
moduli  Ex=50GPa,  E2=8.0GPa,  E2=400GPa  and  Poisson’s  ratios  v,  =  0.30 , 
v2  =  0.28 ,  v3  =  0.20  .  Fet  length  l  =  0.62 m ,  height  h  =  0.05 m  ,  and  crack  length 
a  =  0.0  bn  .  The  left  edge  is  clamped  while  the  right  edge  is  subject  to  compressive  force 
Fx  =  5,0007V  and  constrained  in  the  y-direction.  To  control  for  the  effects  of  bending, 

surface  nodes  at  the  loading  point  are  coupled  in  the  y-direction  thus  preventing  the  free 
end  from  rotating  while  still  allowing  for  axial  displacement. 
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Figure  4.  Composite  beam  with  internal  crack  under  applied  compressive  load. 

2.  FEM  Model  of  Composite  Beam  with  Internal  Crack 

The  same  kind  of  finite  elements  used  in  the  strength  model  are  also  adopted  here 
except  that  in  the  region  of  the  crack  tip  quarter-point  singular  eight-node  elements  (i.e., 
PLANE  183)  are  used  instead.  The  radius  of  the  first  row  elements  is  specified  at 
a! 8  =  1.25 mm  and  the  second  row  ratio  is  set  to  50  percent.  Linear  elastic  plane  stress 
behavior  is  also  assumed  for  this  static  analysis.  A  very  small  opening  (~ a/200) 
facilitates  modeling  of  the  crack  face. 

Each  area  is  freely  meshed  using  smart  size  level  1  elements  with  one  refinement. 
The  right  side  of  the  beam  has  an  applied  total  force  load  Fx  =  -5 kN  distributed  to  each 

node,  constrained  in  the  y  direction.  The  left  side  of  the  beam  is  clamped  and  therefore 
constrained  in  bothx  and  v  directions.  See  Figure  5. 
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Crack  Tip 


Figure  5.  Mesh  of  fracture  mechanics  beam  model  with  internal  crack  (a)  under 
compressive  nodal  force  with  expanded  view  of  crack  tip  ( b ). 


3.  Virtual  Crack  Extension  Method 

In  the  virtual  crack  extension  method,  two  separate  analysis  are  performed,  one 
before  extension  with  crack  length  a  and  the  other  after  extension  with  crack  length 
a  +  Aa. 

In  the  first  analysis,  before  extension,  a  2-D  FEM  model  obtains  the  reaction 
solution  for  resultant  nodal  forces  Fx  and  Fy  at  the  crack  tip.  In  this  case,  the  crack  tip  is 
located  at  node  1412  and  therefore  Fx  and  Fv  are  taken  at  that  node. 

Before  starting  the  second  analysis,  the  crack  length  is  extended.  Extension  of  the 
crack  length  is  carried  out  by  selecting  all  nodes  in  the  vicinity  of  the  crack  tip,  which  is 
defined  here  as  all  nodes  located  within  radial  distance  r  =  a/ 2  =5 mm  of  the  crack  tip, 
and  then  the  selected  nodes  are  scaled  and  moved  in  the  negative  x  direction  for  a 
distance  equivalent  to  1/2  percent  of  the  crack  length.  This  action  results  in  total  crack 
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length  extension  Aa  =  10/;z/;z (0.5%)  =  0.05mm  ,  which  corresponds  to  nodal  scaling 
factor  R  =  40.495/40.5 .  The  crack  tip  remains  at  node  1412,  now  at  an  extended 
distance  from  its  prior  location. 


In  the  second  analysis,  with  crack  length  extension,  a  2-D  FEM  model  obtains  the 
reaction  solution  for  resultant  nodal  displacement  components  ux  and  u  at  the  crack  tip. 

In  this  case,  we  take  the  difference  of  ux  and  uy  between  nearest  neighboring  nodes  1419 
and  1421  post-solution,  using 

(2) 
(3) 


Au.x  Ux-Nl4\9  Ux-Nl42l 


AW,.  Uy-NU19  Uy-N1421 


Linear  interpolation  provides  the  resultant  components  of  virtual  displacement  Aux2  and 
A u  2  at  the  crack  tip,  as  given  by  equations 


AWv2  =  ^~(AUx) 

Ally2=^-{Auy) 


(4) 

(5) 


where  b  is  the  horizontal  distance  between  node  1412  and  either  node  1419  or  node  1421. 
See  Figure  6. 


Figure  6.  Schematic  interpolation  of  crack  tip  forces  &  displacements  for  fracture 
mechanics  beam  model. 
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4.  Energy  Release  Rate 


The  energy  release  rate  G  represents  the  crack’s  driving  force  and  is  the  amount 
of  virtual  work  associated  with  opening  (or  closing)  of  a  crack.  G  is  calculated  for  both 


mode  I  and  mode  II  fractures  using, 


A  a 
A»,2 
A  a 


(6) 

(V) 


where  for  the  constrained  load  at  u  =0  the  crack  surface  is  prevented  from  displacing 


(i.e.,  Am  2  =0)  and  thus  results  in  a  null  energy  release  rate  of  G,  =  0 .  Consequently, 
results  from  the  energy  method  are  presented  here  only  for  mode  II  fractures. 


D.  FRACTURE  MECHANICS  PLATE  MODEL 


The  fracture  mechanics  plate  model  shows  how  an  edge  crack  or  flaw  in  the 
interface  of  a  composite  joint  propagates  under  an  applied  shear  load  and  how  this 
process  is  influenced  by  the  application  of  reinforcing  particles  such  as  CNTs.  These 
analytical  predictions  are  made  by  calculating  fracture  parameters  such  as  energy  release 
rate  in  the  region  of  a  crack  in  order  to  estimate  crack  growth  rate.  While  energy  release 
rate  G  represents  the  amount  of  work  associated  with  a  crack  opening,  it  is  at  the  point  of 
instability  when  the  change  in  potential  energy  equals  the  change  in  the  crack  length  that 
fracture  occurs. 


1.  Composite  Plate  with  Surface  Crack 

Consider  a  composite  plate  with  an  edge  crack  as  shown  in  Figure  7.  The  plate 
comprises  three  linear  elastic  isotropic  materials  having  Young’s  moduli  E]  =50 GPa , 
E2  =  8.0 GPa  ,  E3  =  400 GPa  and  Poisson’s  ratios  v,  =  0.30  ,  v2  =  0.28  ,  v3  =  0.20  .  Let 

length  l  =  0.05m,  height  h=  0.05m,  and  crack  length  a  =  0.01m  .  The  left  edge  is 
clamped  while  the  right  edge  is  constrained  in  the  x  direction  and  subject  to  a  total 
applied  nodal  force  F  =  -5,000N  to  represent  distributed  shearing  load. 
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2.  FEM  Model  of  Composite  Plate  with  Surface  Crack 

A  two-dimensional  FEM  model  in  ANSYS  uses  higher  order  eight-node  elements 
(i.e.,  PLANE82).  Each  node  has  two  degrees  of  freedom,  ux  and  uv,  with  the  degenerated 
element  having  coincident  nodes  K,  L,  and  O  [14]. 

The  crack  tip  uses  quarter-point  singular  eight-node  elements  (i.e.,  PLANE82), 
where  the  radius  of  the  first  row  elements  is  specified  at  a! 8  =  1 25mm  and  the  second 
row  ratio  is  set  to  50  percent.  Linear  elastic  plane  stress  behavior  is  assumed  for  this 
static  analysis.  A  very  small  opening  (~  a/200  )  facilitates  modeling  of  the  crack  face. 

Each  area  is  freely  meshed  using  smart  size  level  1  elements  with  three 
refinements.  The  right  side  of  the  plate  has  applied  total  force  load  Fy  =  -5 kN 

distributed  to  each  node,  constrained  in  the  x  direction.  The  left  side  of  the  plate  is 
clamped  and  therefore  constrained  in  both  x  and  v  directions.  See  Figure  8. 
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(a)  (b) 

Figure  8.  Mesh  of  fracture  mechanics  plate  model  with  surface  crack  (a)  under  shear 
nodal  force  with  expanded  view  of  crack  tip  ( b ). 


3.  Virtual  Crack  Extension  Method 

In  the  virtual  crack  extension  method,  two  separate  analysis  are  performed,  one 
before  extension  with  crack  length  a  and  the  other  after  extension  with  crack  length 
a  +  Aa. 

In  the  first  analysis,  before  extension,  a  2-D  FEM  model  obtains  the  reaction 
solution  for  resultant  nodal  forces  Fx  and  Fy  at  the  crack  tip.  In  this  case,  the  crack  tip  is 
located  at  node  34  and  therefore  Fx  and  Fv  are  taken  at  that  node. 

Before  starting  the  second  analysis,  the  crack  length  is  extended.  Extension  of  the 
crack  length  is  carried  out  by  selecting  all  nodes  in  the  vicinity  of  the  crack  tip,  which  is 
defined  here  as  all  nodes  located  within  radial  distance  r  =  a/2=5mm  of  the  crack  tip, 
and  then  the  selected  nodes  are  scaled  and  moved  in  the  negative  y  direction  for  a 
distance  equivalent  to  1/2  percent  of  the  crack  length.  This  action  results  in  total  crack 
length  extension  Aa  =  10mm (0.5%)  =  0.05mm  ,  which  corresponds  to  nodal  scaling 
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factor  Ry  =3.995/4.000.  Crack  tip  remains  at  node  34,  now  at  an  extended  distance 
from  its  prior  location. 

In  the  second  analysis,  with  crack  length  extension,  a  2-D  FEM  model  obtains  the 
reaction  solution  for  resultant  nodal  displacement  components  ux  and  u  at  the  crack  tip. 

In  this  case,  we  take  the  difference  of  ux  and  u  between  nearest  neighboring  nodes  41 


and  2948  post-solution,  using 


(8) 

(9) 


Linear  interpolation  provides  the  resultant  components  of  virtual  displacement  A ux2  and 


A u  2  at  the  crack  tip,  as  given  by  equations 


A«i2=—  (Ax,)  (10) 

Au,2=-^-(Au,)  (H) 

where  b  is  the  vertical  distance  between  node  34  and  either  of  node  41  or  node  2948.  See 
Figure  9. 
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Figure  9.  Schematic  interpolation  of  crack  tip  forces  &  displacements  for  fracture 
mechanics  plate  model. 


4.  Energy  Release  Rate 


The  energy  release  rate  G  represents  the  crack’s  driving  force  and  is  the  amount 
of  virtual  work  associated  with  opening  (or  closing)  of  a  crack.  G  is  calculated  for  both 
mode  I  and  mode  II  fractures  using, 


1  2  Aa 

1  Aw, 

G'  =  2F’-£ 


(12) 

(13) 


where  for  the  constrained  load  at  ux=  0  the  crack  surface  is  prevented  from  displacing 


(i.e.,  A ux2  =0)  and  thus  results  in  a  null  energy  release  rate  of  G,  =  0 .  Consequently, 
results  from  the  energy  method  are  presented  here  only  for  mode  II  fractures. 


13 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


14 


III.  RESULTS  AND  DISCUSSION 


A.  STRENGTH  MODEL 

The  strength  model  predicts  the  effect  that  reinforcing  particles  have  on  the  state 
of  stress  in  the  interface  of  a  composite  joint.  Several  types  of  models  are  used  for 
comparative  analyses,  including  the  homogeneous  beam,  the  baseline  composite  beam 
without  reinforcing  particles,  and  the  composite  beam  with  reinforcing  particles  where 
different  volume  fractions  are  discretely  tested. 

1.  Homogeneous  Beam 

Starting  with  the  basic  form  of  a  2-D  finite  element  method  (FEM)  model, 
homogeneous  beam  using  base  material  having  Young’s  modulus  E  =  50GPa 
throughout,  produces  nodal  forces  and  displacements  that  result  in  maximum  von  Mises 
equivalent  stress  ae  =  1 5.5MPa  .  See  Table  1. 


Ei 

Ux 

uy 

max 

CT 

max 

Y  max 

e 

max 

V. 

^max 

GPa 

cm 

cm 

kg/cm'S  2 

kg/cm- s  2 

cm/cm 

cm/cm 

kg  -cm  2 /s  2 

50.0 

1.238E-02 

1.520E-04 

31,351 

154,559 

1.63E-04 

3.1 10E-04 

7.857E-01 

Table  1 .  Maximum  stress  and  strain  of  homogeneous  beam  ( E  =  5()GPa  ). 

2.  Baseline  Composite  Beam  without  Reinforcing  Particles 

Next,  using  the  baseline  2-D  FEM  model,  composite  beam  with  base  and  matrix 
material  having  Young’s  moduli  Ex  =  50 GPa  and  E2  =  8.0 GPa  respectively  but  without 
reinforcing  particles,  produces  nodal  forces  and  displacements  that  result  in  maximum 
von  Mises  equivalent  stress  cr,  =  I  1 .5 MPa  in  the  matrix  material.  See  Table  2. 
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Ux 

uy 

^  max 

CT 

max 

Y  max 

e 

max 

V. 

^max 

cm 

cm 

kg/cm-  s  2 

kg/cm- S 

cm/cm 

cm/cm 

kg  ■ cm 2 /s  2 

1.051E-02 

1.587E-03 

17,357 

114,781 

5.550E-04 

1.438E-03 

6.065E+00 

Table  2.  Maximum  stress  and  strain  in  matrix  material  of  baseline  composite  beam 
where  Ei  =  50 GPa  ,  E2  =  8.0 GPa  ,  without  reinforcing  particles. 

The  baseline  model  results  indicate  a  24.4%  decrease  in  the  maximum  von  Mises 
equivalent  stress  <Je  when  compared  to  the  homogeneous  model.  However,  the  opposite 

trend  was  expected  as  the  relatively  softer  matrix  material  now  comprises  greater  than 
eight  percent  by  volume  of  the  composite  beam  thus  lowering  the  overall  beam  stiffness 
creating  a  higher  state  of  internal  stress. 

3.  Composite  Beam  with  Reinforcing  Particles  in  Matrix  Material 

Before  modeling  reinforcing  particles,  three  different  particle  sizes  A  are  used  in 
two  different  volume  fractions  d>  .  Small-,  medium-,  and  large-sized  particles  are  used, 
corresponding  to  particle  size  A  =  1.23mm2 ,  A  =  4.91mm2,  and  A  =  11. 0mm2 
respectively.  The  number  of  particles  varies  in  each  model  in  order  to  hold  constant  the 
volume  fraction  <h  ,  which  is  a  function  of  the  number  of  particles  in  the  matrix  material, 
at  either  of  O  =  7%  or  O  =  14% .  See  Figure  10.  Recall  the  composite  beam  is  clamped 
on  the  left  end  and  has  an  applied  compressive  axial  force  on  the  right. 
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Figure  10.  Composite  beam  matrix  material  containing  reinforcing  particles  of  three 
sizes  A  &  two  volume  fractions  O  :  A0=7%  =  1 .23 mm1  (a),  A<p=l4%  =  1 .23 mm2 

(b),  A0=7%  =4.91mm2  (c),  A^=u%  =  4.9lmm2  (d),  A4l=7%  =  1  I  .Omni1  (e), 
4p=I4%  =11-0 mm2  (f). 


For  the  varying  size  model,  three  distinct  simulations  are  performed,  one  for  each 
particle  size.  The  lowest  number  of  particles  Nd  =  17  and  Nd  =33  use  the  largest  size 

particles  ^4  =  11.0/m/m2  to  maintain  volume  fraction  O  constant.  Similarly,  the  highest 
number  of  particles  Nd  =156  and  Nd  =  296  use  the  smallest  size  A  =  1.23/m/m2  . 

Varying  particle  size  from  A  =  1.23/m/m2  to  A  =  11.0/m/m2  while  holding  constant 
volume  fraction  O  at  either  of  0  =  7%  or  0  =  14%  produces  nodal  forces  and 
displacements  that  result  in  maximum  von  Mises  equivalent  stress  cre  ranging  from 

<je  =  \  62MPa  to  <Je  =  20.1  MPa  in  the  matrix  material.  See  Table  3. 
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Nd 

O 

% 

Ux 

cm 

uy 

cm 

max 

kg/cm- s  2 

a 

max 

kg/cm- s  2 

Y  max 

cm/cm 

e 

max 

cm/cm 

7 / 

^-max 

kg  • cm  2  /s 2 

17 

1.2% 

1.023E-02 

1.275E-03 

31,863 

169,303 

1.020E-03 

2.133E-03 

4.842E+00 

33 

14.0% 

1.012E-02 

1.240E-03 

44,273 

176,270 

1.417E-03 

2.21  IE-03 

3.596E+00 

39 

7.4% 

1.030E-02 

1.386E-03 

42,270 

162,936 

1.353E-03 

2.048E-03 

2.788E+00 

74 

14.0% 

1.014E-02 

1.278E-03 

43,881 

168,642 

1.035E-03 

2.141E-03 

1.875E+00 

156 

7.4% 

1.031E-02 

1.439E-03 

48,700 

181,902 

1.558E-03 

3.224E-03 

9.574E-01 

296 

14.0% 

1.017E-02 

1.309E-03 

48,427 

200,876 

1.550E-03 

2.516E-03 

5.156E-01 

Table  3.  Maximum  stress  and  strain  in  matrix  material  of  composite  beam  where 
Ex=50GPa,  E2=8.0GPa,  E,  =  400 GPa ,  varying  particle  size  from 

A  =  1.23mm1  to  A  =  11. 0mm2  while  holding  constant  volume  fraction  O  at 
either  of  O  =  7%  or  O  =  14%  . 


For  the  varying  particle  size  model,  there  is  an  average  49.3%  increase  in  the 
resultant  ae  for  O  =  7%  and  an  average  58.5%  increase  in  the  resultant  ae  for  O  =14% 
when  compared  to  the  baseline  model.  See  Table  4. 


Particle  Size,  mm 2 

0>=7% 

0=74% 

Small,  1.23 

75.0% 

58.5% 

Medium,  4.91 

42.0% 

46.9% 

Large,  11.0 

47.5% 

53.6% 

Table  4.  Average  increase  in  maximum  vonMises  equivalent  stress  <Je  varying 

reinforcing  particle  size  A  in  matrix  material  of  composite  beam  when 
compared  to  baseline  model. 


Combined  results  show  the  maximum  von  Mises  equivalent  stress  <Je  initially 

decreases  before  rising  slightly  again  as  the  number  of  particles  decreases  with  increasing 
particle  size.  See  Figure  11.  The  reason  for  these  minima  is  not  fully  understood. 
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7%  14% 


Quantity  of  Particles 


Figure  1 1 .  Combined  maximum  von  Mises  equivalent  stresses  ae  in  matrix  material  of 
composite  beam  with  reinforcing  particles  for  volume  fractions  ®  =  7%  & 
®  =  14% . 


The  trend  here  shows  a  higher  volume  fraction  has  a  higher  maximum  von  Mises 
equivalent  stress  cr,  when  compared  to  the  baseline  model,  with  an  overlap  in  stress 

variances  between  the  maximum  and  minimum  in  each  model.  See  Figure  12.  However, 
the  opposite  trend  was  expected  as  the  volume  fraction  goes  up  the  beam  is  generally 
stiffer  and  should  reduce  the  overall  state  of  internal  stress. 
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Figure  12.  Average  maximum  von  Mises  equivalent  stress  ae  in  matrix  material  of 
composite  beam  with  reinforcing  particles  at  volume  fraction  ®  =  7% , 
®  =  14% ;  comparison  to  base  model. 


Overall  these  results  suggest  the  strength  model  is  not  a  reliable  predictor  of  the 
failure  strength  at  the  interface  of  a  composite  beam,  thus  a  different  approach  is  needed. 

B.  FRACTURE  MECHANICS  BEAM  MODEL 

A  fracture  mechanics  based  beam  model  predicts  the  effect  that  reinforcing 
particles  have  on  an  internal  crack  within  the  interface  of  a  composite  joint.  Comparative 
analyses  are  drawn  between  the  baseline  model  having  no  reinforcing  particles  and  a 
composite  beam  with  reinforcing  particles. 

1.  Baseline  Composite  Beam  without  Reinforcing  Particles 

Jumping  right  into  the  baseline  2-D  FEM  model,  composite  beam  with  internal 
crack  but  without  reinforcing  particles,  both  base  and  matrix  material  having  Young’s 
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moduli  El=5QGPa  and  E2=8.0GPa  respectively,  produces  nodal  forces  and 
displacements  at  the  crack  tip  that  result  in  energy  release  rate  Gu  =22.6  mj/m2  .  See 
Table  5. 


Fx 

«x+l 

«x-l 

Aux 

AuX2  Gn 

kg-cm/s  2 

cm 

cm 

cm 

cm  J/m 2 

-65.03 

-9.5869E-03 

-9.5654E-03 

-2.1500E-05 

-3.4677E-06  2.2550E-02 

Table  5.  Energy  release  rate  Gn  in  matrix  material  of  baseline  composite  beam  where 
El  =  50 GPa  ,  E2  =  8.0 GPa  . 

2.  Composite  Beam  with  Reinforcing  Particles  in  Matrix  Material 

Next,  application  of  reinforcing  particles  having  Young’s  modulus  E3  =400 GPa 
to  the  matrix  material  containing  an  internal  crack,  produces  nodal  forces  and 
displacements  at  the  crack  tip  that  result  in  energy  release  rate  Gu  =22.2  mJjm1  .  See 
Table  6. 


Fx 

«x+l 

Ux-l 

Aux 

AUx2  G„ 

kg-cm/s  2 

cm 

cm 

cm 

cm  J/m 2 

-61.73 

-9.4534E-03 

-9.4354E-03 

-1.8000E-05 

-3.6000E-06  2.2222E-02 

Table  6.  Energy  release  rate  Gu  in  matrix  material  of  composite  beam  with  reinforcing 
particles  where  E3  =  400 GPa  . 

Results  from  composite  beam  with  reinforcing  particles  indicate  1.5%  reduction 
in  rate  of  energy  release  Gu  when  compared  to  the  baseline  model.  This  trend  is 
encouraging  and  suggests  the  fracture  mechanics  beam  model  is  a  reliable  predictor  of 
energy  release  rates  for  a  composite  beam  with  internal  crack  at  the  interface. 

Now  that  this  methodology  has  been  validated,  a  more  detailed  fracture 
mechanics  model  is  needed  in  order  to  predict  how  reinforcing  particles  affect  interface 
strength. 
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c. 


FRACTURE  MECHANICS  PLATE  MODEL 


A  fracture  mechanics  based  plate  model  predicts  the  effect  that  reinforcing 
particles  have  on  a  surface  crack  within  the  interface  of  a  composite  joint.  Several  types 
of  models  are  used  for  comparative  analyses,  including  the  homogeneous  plate,  the 
baseline  plate  without  reinforcing  particles,  and  composite  plates  with  reinforcing 
particles  where  the  size,  stiffness,  and  quantity  of  reinforcing  particles  are  discretely 
tested. 


1.  Homogeneous  Plate 

Starting  with  a  basic  form  of  2-D  FEM  model,  the  homogeneous  plate  using  base 
material  having  Young’s  modulus  E  =  50 GP a  throughout,  produces  nodal  forces  and 
displacements  at  the  crack  tip  that  result  in  energy  release  rate  Gu  =8.39  j/m2 .  See 
Table  7. 


Fy 

uy+i 

uy-i 

Ally 

AUy2 

G„ 

kg  ■  cm/s 

cm 

cm 

cm 

cm 

J/m ' 

-3179.4 

-2.0333E-03 

-1.9510E-03 

-8.2300E-05 

-2.6378E-05 

8.3867E+00 

Table  7.  Energy  release  rate  Gn  in  homogeneous  plate  ( E  =  50 GPa  ). 

2.  Baseline  Composite  Plate  without  Reinforcing  Particles 

Next,  the  baseline  2-D  FEM  model,  composite  plate  with  both  base  and  matrix 
material  having  Young’s  moduli  Ex  =50 GPa  and  E2  =8.0 GPa  respectively  but  without 
reinforcing  particles,  produces  nodal  forces  and  displacements  at  the  crack  tip  that  result 
in  energy  release  rate  Gu  =  67.4  J/m2  .  See  Table  8. 


Fy 

uy+i 

uy-i 

Ally 

Ally!  Gjj 

kg-cm/s  2 

cm 

cm 

cm 

cm  J/m 2 

-3606.9 

-4.7695E-03 

-4.1865E-03 

-5.8300E-04 

-1.8686E-04  6.7398E+01 

Table  8.  Energy  release  rate  Gn  in  matrix  material  of  baseline  composite  plate  where 
Ex  =50 GPa  ,  E2  =8.0 GPa  without  reinforcing  particles. 
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The  baseline  model  results  indicate  a  substantial  seven-fold  increase  in  rate  of 
energy  release  Gn  when  compared  to  the  homogeneous  model.  This  trend  is  expected  as 
the  relatively  softer  matrix  material  now  comprises  about  one  third  of  the  composite  plate 
surrounding  the  crack  tip  thus  facilitating  crack  propagation  and  causing  Gn  to  increase 
nearly  one  order  of  magnitude. 

3.  Composite  Plate  with  Reinforcing  Particles  in  Matrix  Material 

Before  modeling  reinforcing  particles,  test  positions  are  established  in  the  matrix 
material  relative  to  the  crack  tip.  Particle  positions  a,  b,  and  c,  are  located  radial  distance 
r  from  the  crack  tip  in  the  negative  x,  negative  y,  and  positive  x  directions  respectively. 
See  Figure  13.  Recall  the  composite  plate  is  clamped  on  the  left  end  and  has  a  transverse 
nodal  shearing  force  on  the  right. 
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Figure  13.  Composite  plate  matrix  material  containing  reinforcing  particles  a,  b,  c  in 
position  radial  distance  r  to  crack  tip. 


4.  Composite  Plate  with  Varying  Particle  Position 

For  the  varying  particle  position  model,  three  distinct  simulations  are  performed, 
one  for  each  particle  position,  followed  by  one  additional  simulation  for  all  three  particles 
together. 
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a.  Particle  at  Position  a 

Varying  particle  a  stiffness  from  E3=8.0GPa  to  E3  =500GPa  while 

holding  constant  particle  size  A  =  0.0396  1/m/m2  at  distance  r  =  1 .  1 24mm  to  crack  tip 
produces  nodal  forces  and  displacements  at  the  crack  tip  that  result  in  energy  release  rates 
ranging  from  Gu  =  67.4  j/m2  down  to  Gu  =  67.3  j/m2 .  See  Table  9. 


E3  Fy  uy+1 

GPa  kg-cm/s 2  cm 

8.0  -3,607.2  -4.7695E-03 

20  -3,606.1  -4.7691E-03 

50  -3,605.2  -4.7689E-03 

80  -3,605.0  -4.7688E-03 

140  -3,604.8  -4.7688E-03 

200  -3,604.7  -4.7688E-03 

250  -3,604.6  -4.7688E-03 

300  -3,604.6  -4.7687E-03 

400  -3,604.5  -4.7687E-03 

500  -3,604.5  -4.7687E-03 


Uy-l  AUy 

cm  cm 

-4.1865E-03  -5.8300E-04 
-4.1864E-03  -5.8270E-04 
-4.1864E-03  -5.8250E-04 
-4.1864E-03  -5.8240E-04 
-4.1864E-03  -5.8240E-04 
-4.1864E-03  -5.8240E-04 
-4.1864E-03  -5.8240E-04 
-4.1864E-03  -5.8230E-04 
-4.1864E-03  -5.8230E-04 
-4.1864E-03  -5.8230E-04 


Auy2  Gn 

cm _ J/m  2 

-1.8686E-04  6.7404E+01 
-1.8676E-04  6.7349E+01 
-1.8670E-04  6.7309E+01 
-1.8667E-04  6.7293E+01 
-1.8667E-04  6.7290E+01 
-1.8667E-04  6.7288E+01 
-1.8667E-04  6.7286E+01 
-1.8663E-04  6.7274E+01 
-1.8663E-04  6.7272E+01 
-1.8663E-04  6.7272E+01 


Table  9.  Energy  release  rate  Gn  for  reinforcing  particle  a  in  matrix  material  of 
composite  plate  varying  particle  stiffness  from  E 3  =  8.0 GPa  to  E3  =  500 GPa 

with  particle  size  A  =  0.03961mm2  at  distance  r  =  1.124mm  to  crack  tip. 

b.  Particle  at  Position  b 

Varying  particle  b  stiffness  from  E3=8.0GPa  to  E3  =500GPa  while 

holding  constant  particle  size  A  =  0.03961mm2  at  distance  r  =  1.124/m/m  to  crack  tip 
produces  nodal  forces  and  displacements  at  the  crack  tip  that  result  in  energy  release  rates 
ranging  from  Gu  =  67 Aj/m2  down  to  Gu  =  66.8  j/m1 .  See  Table  10. 
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E3 

GPa 

Fy  Uy+1 

kg-cm/s 2  cm 

Uy_, 

cm 

AUy 

cm 

AUy2 

cm 

G„ 

J/m  2 

8.0 

-3,606.3  -4.7694E-03 

-4.1865E-03 

-5.8290E-04 

-1.8683E-04 

6.7375E+01 

20 

-3,598.8  -4.7685E-03 

-4.1868E-03 

-5.8170E-04 

-1.8644E-04 

6.7097E+01 

50 

-3,594.5  -4.7680E-03 

-4.1869E-03 

-5.81 10E-04 

-1.8625E-04 

6.6948E+01 

80 

-3,593.2  -4.7678E-03 

-4.1869E-03 

-5.8090E-04 

-1.8619E-04 

6.6900E+01 

140 

-3,592.1  -4.7677E-03 

-4.1870E-03 

-5.8070E-04 

-1.8612E-04 

6.6857E+01 

200 

-3,591.7  -4.7677E-03 

-4.1870E-03 

-5.8070E-04 

-1.8612E-04 

6.6849E+01 

250 

-3,591.5  -4.7676E-03 

-4.1870E-03 

-5.8060E-04 

-1.8609E-04 

6.6834E+01 

300 

-3,591.3  -4.7676E-03 

-4.1870E-03 

-5.8060E-04 

-1.8609E-04 

6.6830E+01 

400 

-3,591.1  -4.7676E-03 

-4.1870E-03 

-5.8060E-04 

-1.8609E-04 

6.6827E+01 

500 

-3,591.0  -4.7676E-03 

-4.1870E-03 

-5.8060E-04 

-1.8609E-04 

6.6825E+01 

Table  10.  Energy  release  rate  Gn  for  reinforcing  particle  b  in  matrix  material  of 
composite  plate  varying  particle  stiffness  from  Zs3  =  8.0 GPa  to  E,  =  500 GPa 

with  particle  size  A  =  0.03961mm2  at  distance  r  =  1.124mm  to  crack  tip. 


c.  Particle  at  Position  c 

Varying  particle  c  stiffness  from  Ei=8.0GPa  to  E3  =  500GPa  while 

holding  constant  particle  size  A  =  0.0396 1mm2  at  distance  r  =  1.124mm  to  crack  tip 
produces  nodal  forces  and  displacements  at  the  crack  tip  that  result  in  energy  release  rates 
ranging  from  Gn  =  67.4  j/m2  down  to  Gu  =  67.2  j/m2  .  See  Table  11. 
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E3 

GPa 

Fy  Uy+1 

kg-cm/s 2  cm 

Uy_, 

cm 

AUy 

cm 

AUy2 

cm 

G„ 

J/m  2 

8.0 

-3,607.1  -4.7695E-03 

-4.1865E-03 

-5.8300E-04 

-1.8686E-04 

6.7402E+01 

20 

-3,603.4  -4.7693E-03 

-4.1866E-03 

-5.8270E-04 

-1.8676E-04 

6.7298E+01 

50 

-3,601.0  -4.7691E-03 

-4.1866E-03 

-5.8250E-04 

-1.8670E-04 

6.7230E+01 

80 

-3,600.3  -4.7691E-03 

-4.1867E-03 

-5.8240E-04 

-1.8667E-04 

6.7206E+01 

140 

-3,599.8  -4.7690E-03 

-4.1867E-03 

-5.8230E-04 

-1.8663E-04 

6.7185E+01 

200 

-3,599.5  -4.7690E-03 

-4.1867E-03 

-5.8230E-04 

-1.8663E-04 

6.7179E+01 

250 

-3,599.4  -4.7690E-03 

-4.1867E-03 

-5.8230E-04 

-1.8663E-04 

6.7177E+01 

300 

-3,599.3  -4.7690E-03 

-4.1867E-03 

-5.8230E-04 

-1.8663E-04 

6.7175E+01 

400 

-3,599.2  -4.7690E-03 

-4.1867E-03 

-5.8230E-04 

-1.8663E-04 

6.7174E+01 

500 

-3,599.2  -4.7690E-03 

-4.1867E-03 

-5.8230E-04 

-1.8663E-04 

6.7174E+01 

Table  11.  Energy  release  rate  G//  for  reinforcing  particle  c  in  matrix  material  of 
composite  plate  varying  particle  stiffness  from  E 3  =  8.0 GPa  to  E,  =  500 GPa 

with  particle  size  A  =  0.03961mm2  at  distance  r  =  1.124mm  to  crack  tip. 


d.  Particles  at  Positions  a,  b,  c 

Varying  combined  particles  abc  stiffness  from  E:,  =  S.OGPa  to 
Ei=500GPa  while  holding  constant  particle  size  A  =  0.0396 1mm2  at  distance 
r  =  1.124 mm  to  crack  tip  produces  nodal  forces  and  displacements  at  the  crack  tip  that 
result  in  energy  release  rates  ranging  from  G„  =  67.4  j/ nr  down  to  Gu  =  66.5  j/m2 . 
See  Table  12. 
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E3 

GPa 

Fy 

kg-cm/s  2 

Uy+i 

cm 

Uy-l 

cm 

AUy 

cm 

Auy2 

cm 

G„ 

J/m  2 

8.0 

-3606.8 

-4.7695E-03 

-4.1865E-03 

-5.8300E-04 

-1.8686E-04 

6.7396E+01 

20 

-3594.3 

-4.7679E-03 

-4.1868E-03 

-5.81 10E-04 

-1.8625E-04 

6.6944E+01 

50 

-3586.8 

-4.7670E-03 

-4.1870E-03 

-5.8000E-04 

-1.8590E-04 

6.6678E+01 

80 

-3584.5 

-4.7668E-03 

-4.1870E-03 

-5.7980E-04 

-1.8583E-04 

6.6612E+01 

140 

-3582.6 

-4.7665E-03 

-4.1871E-03 

-5.7940E-04 

-1.8571E-04 

6.6531E+01 

200 

-3581.8 

-4.7664E-03 

-4.1871E-03 

-5.7930E-04 

-1.8567E-04 

6.6504E+01 

250 

-3581.4 

-4.7664E-03 

-4.1871E-03 

-5.7930E-04 

-1.8567E-04 

6.6497E+01 

300 

-3581.2 

-4.7664E-03 

-4.1871E-03 

-5.7930E-04 

-1.8567E-04 

6.6493E+01 

400 

-3580.9 

-4.7663E-03 

-4.1871E-03 

-5.7920E-04 

-1.8564E-04 

6.6476E+01 

500 

-3580.6 

-4.7663E-03 

-4.1871E-03 

-5.7920E-04 

-1.8564E-04 

6.6471E+01 

Table  12.  Energy  release  rate  Gn  for  reinforcing  particle  abc  in  matrix  material  of 
composite  plate  varying  particle  stiffness  from  E3  =  8.0 GPa  to  E:,  =  500 GPa 

with  particle  size  A  =  0.0396 Imm2  at  distance  r  =  1.124mm  to  crack  tip. 


For  the  discrete  varying  particle  stiffness  models,  there  is  an  average  0.4% 
decrease  in  the  resultant  energy  release  rate  Gn  when  compared  to  the  baseline  model. 
This  average  increases  to  1.2%  when  all  three  positions  are  combined.  See  Table  13. 
Individual  results  at  E}=8.0GPa  for  each  varying  stiffness  model  demonstrate 
continuity  with  the  baseline  model. 


Particle  Location 

Gn 

a 

-0.2% 

b 

-0.8% 

c 

-0.3% 

abc 

-1.2% 

Table  13.  Average  decrease  in  energy  release  rate  Gn  for  reinforcing  particles  a,b,c, 
when  compared  to  baseline  model. 


Combined  results  show  energy  release  rate  Gn  decreases  as  particle 
stiffness  increases  when  compared  to  the  baseline  model.  See  Figure  14.  In  addition,  Gn 
asymptotically  tends  to  zero  for  elastic  modulus  above  E3  =  200GPa 
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Modulus  of  Elasticity  E3  (GPa) 


Figure  14.  Combined  energy  release  rates  Gn  for  reinforcing  particles  a,  b,  c,  in  matrix 
material  of  composite  plate  varying  particle  stiffness  from  E3  =  8.0 GPa  to 

E3  =  500GPa  with  particle  size  A  =  0.0396 Imm2  at  distance  r  =  1.124mm 
to  crack  tip. 


The  trend  here  shows  that  while  particle  b  has  the  largest  single  decrease 
in  Gn,  a  more  favorable  decrease  is  predicted  when  all  three  positions  a,  b,  and  c,  are 
combined.  See  Figure  15.  These  results  are  expected  as  particle  b  is  positioned 
favorably  to  have  greatest  influence  over  the  amount  of  energy  released  in  mode  II 
fractures  of  this  kind. 
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Figure  15.  Average  energy  release  rate  Gn  for  reinforcing  particles  a,  b,  c,  in  matrix 
material  varying  particle  stiffness  from  E3  =8.0 GPa  to  E3  =500 GPa  with 

particle  size  A  =  0.0396 Imm2  at  distance  r  =  1.124mm  to  crack  tip; 
comparison  to  baseline  model. 

5.  Composite  Plate  with  Varying  Particle  b  Distance  to  Crack  Tip 

Building  on  prior  results,  particle  b  is  again  selected  in  varying  distances  to  the 
crack  tip.  In  this  case,  varying  particle  b  distance  r  to  crack  tip  from  r  =  0  to 
r  =  2.937mm  while  holding  constant  particle  stiffness  E3=400GPa  produces  nodal 
forces  and  displacements  at  the  crack  tip  that  result  in  energy  release  rates  ranging  from 
Gu  =  62.5  jjm1  up  to  Gn  =  67.3  j/m2 .  See  Table  14. 
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r 

cm 

Fy 

kg-cm/s  2 

uy+1 

cm 

Uy.i 

cm 

AUy 

cm 

AUy2 

cm 

G„ 

J/m  2 

0.0000 

-3481.9 

-4.7556E-03 

-4.1957E-03 

-5.5990E-04 

-1.7946E-04 

6.2484E+01 

0.0312 

-3537.5 

-4.7630E-03 

-4.1902E-03 

-5.7280E-04 

-1.8359E-04 

6.4945E+01 

0.0625 

-3574.7 

-4.7661E-03 

-4.1879E-03 

-5.7820E-04 

-1.8532E-04 

6.6247E+01 

0.0938 

-3591.2 

-4.7676E-03 

-4.1870E-03 

-5.8060E-04 

-1.8609E-04 

6.6829E+01 

0.1124 

-3591.1 

-4.7676E-03 

-4.1870E-03 

-5.8060E-04 

-1.8609E-04 

6.6827E+01 

0.1428 

-3592.1 

-4.7675E-03 

-4.1868E-03 

-5.8070E-04 

-1.8612E-04 

6.6857E+01 

0.1625 

-3597.4 

-4.7682E-03 

-4.1867E-03 

-5.8150E-04 

-1.8638E-04 

6.7048E+01 

0.1857 

-3602.1 

-4.7686E-03 

-4.1864E-03 

-5.8220E-04 

-1.8660E-04 

6.7216E+01 

0.2383 

-3603.5 

-4.7687E-03 

-4.1863E-03 

-5.8240E-04 

-1.8667E-04 

6.7265E+01 

0.2937 

-3603.6 

-4.7687E-03 

-4.1862E-03 

-5.8250E-04 

-1.8670E-04 

6.7279E+01 

Table  14.  Energy  release  rate  Gu  for  reinforcing  particle  b  in  matrix  material  of 
composite  plate  varying  distance  r  to  crack  tip  from  r  =  0  to  r  =  2.931mm 
with  particle  stiffness  E3  =  400 G Pa  . 


For  the  varying  particle  b  distance  model,  there  is  an  average  1.6%  decrease  in  the 
resultant  energy  release  rate  Gn  when  compared  to  the  prior  baseline  models.  See  Table 
15.  In  addition,  individual  results  from  r  =  1.124mm  demonstrate  continuity  with  the 
prior  stiffness  model. 


Distance  to 
Crack  Tip,  mm 

Gn 

0 

-7.3% 

0.3 

-3.6% 

0.6 

-1.7% 

0.9 

-0.9% 

1.1 

-0.9% 

1.4 

-0.8% 

1.6 

-0.5% 

1.9 

-0.3% 

2.4 

-0.2% 

2.9 

-0.2% 

Table  15.  Average  decrease  in  energy  release  rate  Gn  for  particle  b  varying  distance  r  to 
crack  tip  when  compared  to  baseline  model. 
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Combined  results  show  energy  release  rate  Gu  decreases  as  particle  b  distance  to 
the  crack  tip  decreases  (i.e.,  closest  to  the  crack  tip).  See  Figure  16.  In  addition,  Gu 
asymptotically  tends  to  zero  at  distances  greater  than  r  =  2mm  ,  providing  a  clear  line  of 
demarcation  where  reinforcing  particles  no  longer  influence  crack  propagation. 


Figure  16.  Energy  release  rate  Gu  for  reinforcing  particle  b  in  matrix  material  of 
composite  plate  varying  distance  r  to  crack  tip  from  r  =  0  to  r  =  3 mm  with 
particle  stiffness  E3  =  400 GPa  . 

6.  Composite  Plate  with  Varying  Particle  Size  at  Prescribed  Distance  to 
Crack  Tip 

Particle  b  is  again  used  for  varying  particle  size  at  prescribed  distance  r  to  crack 
tip.  In  this  case,  varying  particle  b  size  from  A  =  0.01  mm2  to  A  =  2 mm2  while  holding 
constant  particle  stiffness  E3  =  400C//C/  at  r  =  \.\ 24mm  to  crack  tip  produces  nodal 
forces  and  displacements  at  the  crack  tip  that  result  in  energy  release  rates  ranging  from 
Gu  =  67.2  j/m1  down  to  Gn  =  53.6  j/m2 .  See  Table  16. 
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A 

2 

mm 

Fy 

kg-cm/s 

Uy+1 

cm 

Uy-l 

cm 

AUy 

cm 

Auy2 

cm 

G„ 

J/m  ' 

0.01 

-3602.2 

-4.7689E-03 

-4.1866E-03 

-5.8230E-04 

-1.8663E-04 

6.7230E+01 

0.04 

-3591.1 

-4.7676E-03 

-4.1870E-03 

-5.8060E-04 

-1.8609E-04 

6.6827E+01 

0.15 

-3548.8 

-4.7625E-03 

-4.1884E-03 

-5.7410E-04 

-1.8401E-04 

6.5300E+01 

0.36 

-3481.6 

-4.7538E-03 

-4.1901E-03 

-5.6370E-04 

-1.8067E-04 

6.2903E+01 

0.70 

-3390.1 

-4.7409E-03 

-4.1914E-03 

-5.4950E-04 

-1.7612E-04 

5.9707E+01 

1.17 

-3298.4 

-4.7264E-03 

-4.1914E-03 

-5.3500E-04 

-1.7147E-04 

5.6559E+01 

1.82 

-3210.7 

-4.7103E-03 

-4.1893E-03 

-5.2100E-04 

-1.6699E-04 

5.3615E+01 

Table  16.  Energy  release  rate  Gn  for  reinforcing  particle  b  in  matrix  material  of 
composite  plate  varying  particle  size  from  A  =  0.01mm2  to  A  =  2mm1  with 
particle  stiffness  Z2  =  AOOGPa  at  prescribed  distance  r  =  1.124mm  to  crack 
tip. 


For  the  varying  particle  size  model,  there  is  an  average  8.4%  decrease  in  the 
resultant  energy  release  rate  Gn  when  compared  to  the  prior  baseline  model.  See  Table 
17.  In  addition,  individual  results  from  A  =  0.04mm2  demonstrate  continuity  with  the 
prior  distance  model. 


Particle  Size,  mm' 

Gn 

0.01 

-0.3% 

0.04 

-0.8% 

0.15 

-3.1% 

0.36 

-6.7% 

0.70 

-11% 

1.2 

-16% 

1.8 

-21% 

Table  17.  Average  decrease  in  energy  release  rate  Gn  for  particle  b  while  varying  size  A 
when  compared  to  baseline  model. 


Combined  results  show  energy  release  rate  Gn  monotonically  decreases  as  particle 
b  size  increases.  See  Figure  17. 
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Figure  17.  Energy  release  rate  Gu  for  reinforcing  particle  b  in  matrix  material  of 
composite  plate  varying  particle  size  from  A  =  0.01mm2  to  A  =  2mm1  with 
particle  stiffness  E 3  =400 GPa  at  prescribed  distance  r  =  1 . 124mm  to  crack 
tip. 

7.  Composite  Plate  with  Varying  Particle  Stiffness  Nearest  Crack  Tip 

Next,  particle  b  is  located  near  the  crack  tip  at  distance  r  =  0.3 mm  and  under 
varying  particle  stiffness  from  E 3  =  8.0 GPa  to  E3  =  500 GPa  produces  nodal  forces  and 
displacements  at  the  crack  tip  that  result  in  energy  release  rates  ranging  from 
Gu  =67 .3  j/m1  down  to  Gu  =64.9  j/m2  .  See  Table  18. 
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E3 

GPa 

Fy  Uy+1 

kg-cm/s 2  cm 

Uy-l 

cm 

AUy 

cm 

AUy2 

cm 

G„ 

J/m  2 

8.0 

.3604.0  -4.7692E-03 

-4.1867E-03 

-5.8250E-04 

-1.8670E-04 

6.7286E+01 

20 

-3571.3  -4.7662E-03 

-4.1884E-03 

-5.7780E-04 

-1.8519E-04 

6.6138E+01 

50 

-3553.0  -4.7645E-03 

-4.1894E-03 

-5.7510E-04 

-1.8433E-04 

6.5491E+01 

80 

-3547.3  -4.7640E-03 

-4.1897E-03 

-5.7430E-04 

-1.8407E-04 

6.5295E+01 

140 

-3542.5  -4.7635E-03 

-4.1899E-03 

-5.7360E-04 

-1.8385E-04 

6.5127E+01 

200 

-3540.4  -4.7633E-03 

-4.1900E-03 

-5.7330E-04 

-1.8375E-04 

6.5055E+01 

250 

-3539.3  -4.7632E-03 

-4.1901E-03 

-5.7310E-04 

-1.8369E-04 

6.5012E+01 

300 

-3538.5  -4.763 IE-03 

-4.1901E-03 

-5.7300E-04 

-1.8365E-04 

6.4986E+01 

400 

-3537.5  -4.7630E-03 

-4.1902E-03 

-5.7280E-04 

-1.8359E-04 

6.4945E+01 

500 

-3536.9  -4.7629E-03 

-4.1902E-03 

-5.7270E-04 

-1.8356E-04 

6.4923E+01 

Table  18.  Energy  release  rate  Gn  for  reinforcing  particle  b  in  matrix  material  of 
composite  plate  varying  particle  stiffness  from  E 3  =  8.0 GPa  to  E3  =  500 GPa 
with  particle  located  nearest  the  crack  tip  at  distance  r  =  0.3 mm  . 


Results  from  varying  particle  stiffness  nearest  the  crack  tip  show  there  is  an 
average  2.9%  decrease  in  the  resultant  energy  release  rate  Gn  when  compared  to  the  prior 
baseline  model.  See  Table  19.  In  addition,  individual  results  from  E3  =S.0GPa 
demonstrates  continuity  with  the  baseline  model. 


Particle  Stiffness 

Gn 

1.0 

-0.2% 

2.5 

-1.9% 

6.3 

-2.8% 

10.0 

-3.1% 

17.5 

-3.4% 

25.0 

-3.5% 

31.3 

-3.5% 

37.5 

-3.6% 

50.0 

-3.6% 

62.5 

-3.7% 

Table  19.  Average  decrease  in  energy  release  rate  Gn  for  particle  b  while  varying 
stiffness  E-,  when  compared  to  baseline  model. 
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Combined  results  show  energy  release  rate  Gu  decreases  as  particle  b  stiffness 
increases.  See  Figure  18. 


Figure  18.  Energy  release  rate  Gu  for  reinforcing  particle  b  in  matrix  material  of 
composite  plate  varying  particle  stiffness  from  E2  =  8.0  GPa  to 

E3  =  500 GPa  with  particle  located  nearest  the  crack  tip  at  r  =  0.3 mm  . 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

Two-dimensional  continuum  modeling  reliably  and  efficiently  predict  how  the 
size,  quantity,  and  stiffness  of  reinforcing  particles  such  as  carbon  nanotubes  (CNTs) 
affect  failure  mechanisms  at  the  interface  of  composite  structures,  depending  on  which 
methodology  is  used. 

The  strength  model  revealed  that  adding  reinforcing  particles  to  the  matrix 
material  of  a  scarf  joint  causes  the  maximum  von  Mises  equivalent  stress  to  increase. 
This  result  was  contrary  to  what  had  been  expected.  Adding  reinforcing  particles  to  the 
interface  of  the  composites  should  have  increased  the  interface  strength.  Therefore, 
based  on  these  results,  the  strength  model  does  not  reliably  predict  the  effect  that 
reinforcing  particles  have  on  state  of  stress  in  the  interface  of  a  composite  joint. 

The  fracture  mechanics  based  beam  model  illustrated  that  application  of 
reinforcing  particles  to  the  matrix  material  of  a  composite  beam  with  an  internal  crack 
causes  the  energy  release  rate  to  favorably  decrease  approximately  two  percent.  This 
result  was  encouraging  and  served  as  the  basis  for  validation  of  the  methodology  to 
proceed  with  a  more  detailed  fracture  mechanics  model  of  the  interface. 

The  fracture  mechanics  plate  model  revealed  many  aspects  of  particle 
reinforcement  in  a  composite  joint.  First,  application  of  three  reinforcing  particles  near 
the  crack  tip  favorably  reduces  the  energy  release  rate  one  percent  on  average.  When 
presented  discretely  the  magnitudes  of  one  particle  were  greater  than  the  others,  thus 
energy  release  rate  is  a  strong  function  of  particle  position  relative  to  the  crack  tip. 
Second,  a  single  reinforcing  particle  located  at  the  crack  tip  favorably  reduces  the  energy 
release  rate  seven  percent  on  average,  thus  energy  release  rate  is  a  strong  function  of 
reinforcing  particle  distance  to  the  crack  tip.  Third,  energy  release  rate  monotonically 
decreases  as  the  reinforcing  particle  size  increases.  Finally,  application  of  reinforcing 
particles  that  are  ten  times  stiffer  than  the  surrounding  matrix  material  favorably  reduces 
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the  energy  release  rate  by  3  percent  on  average.  Therefore,  based  on  these  combined 
results,  the  fracture  mechanics  plate  model  appears  to  reliably  predict  the  effect  that 
reinforcing  particles  have  on  energy  release  rate  in  the  interface  of  a  composite  plate  with 
surface  crack.  Additionally,  model  setup  time  is  estimated  to  be  approximately  one  hour 
with  solution  time  of  less  than  one  minute  for  each  iteration.  Moreover,  model  storage 
requirements  are  nominal  with  file  size  of  approximately  12MB,  which  includes  all  pre¬ 
processing  and  meshing.  Therefore,  in  addition  to  being  reliable,  the  model  is  also  highly 
efficient  in  terms  of  time  and  space. 

B.  RECOMMENDATIONS 

Continuum  modeling  of  interface  failure  represents  reliable  and  efficient 
modeling  of  composite  interfaces  designed  to  improve  their  failure  strength  through  the 
addition  of  nanoscale  reinforcing  particles  such  as  carbon  nanotubes  (CNTs).  This  type 
of  modeling  is  a  useful  tool  to  predict  and  improve  the  failure  strength  of  composite 
structures.  The  model  can  be  used  to  help  focus  future  research  in  the  structural 
applications  of  CNTs,  especially  within  testing  and  evaluation  of  CNTs  in  composite 
scarf  joint  interfaces.  In  addition,  atomic  level  modeling  is  recommended  for  the 
interface  failure  to  be  compared  with  the  continuum  models.  The  atomic  model  can 
include  nanoparticles  such  as  CNTs  as  well  as  the  polymer  matrix  as  they  are  without 
smearing.  Additionally,  multi-scale  modeling  will  be  beneficial  to  include  different 
length-scale  characteristics  in  the  design  and  analysis  of  composite  joints  containing 
nanoparticles. 
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